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Introduction  and  Summary 

During  the  first  quarter  of  this  contract,  seismic  source 
depth  determination  techniques  (developed  during  the  two  previous 
contracts),  were  formulated  into  an  automated  analysis  procedure. 

In  this  procedure,  the  source  depth  is  determined  by  the  degree 
to  which  cepstrum  patterns,  computed  from  different  portions  of 
multi - station  data,  agree  with  those  expected  for  a given  source 
depth.  The  variation  in  differential  travel  times  caused  by  the 
location  of  stations  and  the  presence  of  later  propagation  modes 
are  both  accounted  for,  and  the  depth  phase  information  in  the 
entire  seismogram  is  involved  in  the  depth  estimate.  The  analysis 
procedure  involves  the  computation  of  the  cepstrum,  cepstrum 
match  filter  outputs,  stochastic  averaging,  and  use  of  the 
differential  travel  time  for  several  propagation  modes.  The 
travel  time  data,  automatically  accessed  in  the  analysis  procedure, 
is  stored  as  coefficients  defining  a polynomial  surface  fitted  to 
the  three-dimensional  data  of  travel  time  versus  source  depth  and 
source  to  receiver  epicenter  distance. 

In  the  current  procedure,  multi-channel  seismic  data  is  read 
in  and  a plot  of  the  cumulative  CMF  outputs  versus  source  depth  is 
computed.  The  Illinois  Event  of  11/9/68  has  been  analyzed  thus 
far,  and  the  output  indicates  a clear  detection  of  a 25  km  source 
depth,  correct  for  this  event.  We  now  discuss  in  more  detail  the 
analysis  procedure,  the  computer  storage  of  the  travel  time 
information,  and  the  application  of  the  seismic  source  depth 
analysis . 


Major  Accoir.i'lishments 
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1 . Seismic  Source  Depth* Analysis  Procedure 

The  seismic  source  depth  analysis  procedure  is  illustrated 
in  the  flow  chart  of  Figure  1.  The  seismic  data  recorded  over  a 
suite  of  stations  are  selected  to  start  at  the  onset  of  the  P wave. 
Next  the  data  sample  length,  governing  the  maximum  differential 
delay  time  observable,  and  the  length  of  coda  to  be  analyzed, 
are  selected. 

At  this  point,  a set  of  consecutive  data  samples  covering  the 
selected  coda  length  for  each  station  recording  has  been  deter- 
mined. A data  sample  is  selected  and  a cepstrum  and  cepstrum 
matched  filter  output  is  calculated.  The  procedure  used  to  calcu- 
late che  CMF  is  illustrated  in  Figure  2.  The  first  of  a range  of 
trial  $ource  depths  a^-''  then  selected  for  the  analysis.  For  this 
trial  depth  and  station  epicenter  distance,  the  differential 
travel  times  for  the  following  propagation  modes  are  accessed 
from  storage:  pP-P,  PP-P,  pPP-PP,  PPP-P,  pPPP-PPP,  PcP-P, 

pPcP-PcP.  For  tiie  start  and  end  time  of  the  data  sample,  it  is 
decided  which  modes  will  contribute  to  the  cepstrum  for  this 
particular  data  sample.  The  maximum  CMF  output  within  a given 
time  window  of  the  expected  time  delay  between  the  surface 
reflected  arrival  and  the  direct  arrival  for  each  of  these  pro- 
pagation modes  is  accumulated  in  the  output  amplitude  for  this 
trial  source  depth.  For  this  same  sample,  the  procedure  is 
repeated  over  the  range  of  trial  depths  used.  This  completes 
the  analysis  of  this  data  sample  and  the  next  data  sample  is 
likewise  processed.  One  then  plots  the  cumulative  CMF  output 
as  a function  of  depth,  and  a peak  in  the  distribution  will 
indicate  the  depth  estimate. 
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FIGURE  1.  SF.ISKIC  SOURCE  DEPTH  DUTI! RUINATION  PROCEDURE 
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FIGURE  2.  CALCULATION  OF  CEPSTRUM  MATCHED  FILTER  (CMF)  OUTPUT 
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This  entire  procedure  is  then  repeated  for  different  sample 
lengths  and/or  coda  lengths. 

The  procedure  was  applied  to  the  Illinois  Event  (11/9/68) 
using  a sample  length  of  256  points  and  five  consecutive  coda 
samples  for  each  of  the  six  station  recordings.  The  resulting 
output  as  a function  of  depth  (in  km)  is  plotted  in  Figure  3 and 
a clear  detection  of  the  25  km  source  depth  for  the  event  is 
evident . 

In  the  next  section  we  discuss  the  techniques  developed  to 
facilitate  the  storage  and  accessing  of  the  differential  travel 
times  used. 


Polynomial  Surface  Representation  of  Differential  Travel  Times 


To  account  for  depth  phase  delay  time  variations  along  seismic 
coda,  one  needs  access  to  the  differential  travel  times  for 
various  seismic  phases.  From  previous  work,  these  travel  time 
differences  were  obtained  by  the  application  of  a ray  tracing 
program  based  on  the  spherically  symmetric  isotropic  earth  velocity 
model  used  for  the  BSSA  seismological  tables.  For  this  current 
work,  an  analytic  three-dimensional  representation  of  these  data 
was  developed  in  order  to  facilitate  computer  access  and  to  per- 
form the  necessary  interpolation  of  these  values.  We  now  describe 
the  procedure  for  obtaining  the  representation  of  these  differential 
travel  times  as  a function  of  source  depth  and  source  to  receiver 
distance  for  the  following  seismic  propagation  modes:  pP-P, 

PP-P,  pPP-FP,  PPP-P,  pPPP-PPP,  PeP-P,  pPeP-PeP. 
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HICURG  3.  KHSULTS  OP  SOURCI:  DEPTH  ANALYSIS  PO'>  ILLINOIS  Evn.i- 
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Case  I:  Surface  Fits  to  pP-P,  pPPP-PPP,  pPP-pp,  pPcP-PcP 


An  examination  of  Figure  4 suggests  we  represent  these  data 
as  a double  power  series  in  which  the  travel  time  difference  t 
is  the  dependent  variable  and  the  depth  d and  epicenter  distances 
A are  the  independent  variables.  The  surface  of  travel  time 
differences  is  then  given  by 

(d.A)  * ^ T.  . A^d^ 

^ ij 
11 

The  coefficients  and  the  number  of  them  needed  were  determined 
as  follows.  Each  curve  at  constant  A^  was  found  to  be  adequately 
represented  by  a cubic  in  d; 

T(d,Ac)  “ 

Here  t vanishes  at  zero  depth.  The  least  squares  values  of  t.(a) 
were  then  in  turn  adjusted  to  a power  series  in  L\ 

N 

T'iCA)  - ^ T , i = l,2,3 

j=0  ^ 

Reliable  input  data  for  a less  than  10°  were  not  available,*  so  no 
attempt  was  made  to  force  t.(a)  to  vanish  at  zero  epicenter  distance. 
Thus,  our  representation  should  be  used  only  for  L>,  10°.  All  the 
data  sets  of  Case  I can  be  represented  by  N=9.  Examples  of  the 
effectiveness  of  this  representation  are  shown  in  Figures  5a 
through  5d  for  pP-P,  pPP-PP,  pPPP-ppp,  and  pPeP-PeP. 
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Plot  of  travel  time  differ 
ence  (seconds)  for  pP-P 
vs.  depth  (km)  for  given 
epicenter  distance.'  (A®) 


Depth  (km) 
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Each  differential  travel  time  surface  considered  in  this 
work,  except  those  representing  core  reflections,  exhibit  an 
unallowed  region,  i.e.,  there  exists  a range  of  epicenter  distances 
for  which  below  a certain  depth  the  surface  reflected  mode  is  not 
received.  The  boundary  of  such  a region  projected  intc  the  (d  a) 
plane  for  pP-P  is  shown  in  Figure  6. 


FIGURE  6 


These  boundaries  can  be  adequately  represented  by  a parabola  so 
that  our  polynomial  model  is  subje:t  to  the  constraint  that  for 
a given  A we  must  have  d < a+bA^,  where  a and  b have  been  deter- 
mined from  least  square  adjustments  to  plots  like  that  of  Figure  6. 

In  Figure  5,  the  unallowed  region  is  shown  as  asterisk  table 
entries. 


Case  II:  Line  Fits  to  PP-P,  ppp-p,  PcP-p 


These  travel  time  difference  plots  ca  - be  represented  as  s ngle 
curves  depending  on  a only,  since  between  5 and  100  km,  the  depth 
depenilence  is  very  weak.  Thus  we  have 
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T(d,A)  - T (a)  «y  T .A^. 

0 Zj  OJ 

j-0 


The  PP-P  and  PPP-P  differences  appear  to  vanish  as  A-^o,  so  for 
these  survey;  we  set  - 0 and  also  find  N«=8  is  sufficient.  For 
FcP-P,  the  core  reflection  insures  that  t *f  0.  Here  N»6  is 
sufficient.  °° 

This  polynomial  surface  representation  of  the  differential 
travel  times  can  be  stored  using  a tot’l  of  143  coefficients  and 
requires  less  than  .02  seconds  to  compute  these  times  for  eight 
propagation  modes  on  the  CDC  6600  computer. 


Future  Plans 


A prototype  of  the  automated  seismic  source  depth  analysis 
procedure  has  been  developed  and  applied  to  the  first  minute  of 
the  seismograms  recorded  for  the  Illinois  Event.  During  the  next 
quarter  we  will  apply  this  analysis  to  succeeding  data  in  which  the 
propagation  modes  PP,ppp  and  PcP  are  present.  We  will  investigate 
the  effectiveness  of  this  analysis  in  enabling  constructive  use 
of  the  depth  phase  information  and  to  determine  what  modifications 
may  be  necessary.  We  will  also  investigate,  through  the  application 
of  this  analysis  to  additional  events,  those  modifications  nec- 
essary to  enhance  the  statistical  significance  and  resolution  of 
these  source  depth  estimates.  During  this  period,  we  will  also 
begin  to  investigate  the  usefulness  of  the  Maximum  Entropy  tech- 
niques in  extending  the  applicability  of  the  analysis  procedures 
to  shallow  events. 
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